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ABSTRACT 

We revisit the issue of the constancy of the dark matter (DM) and baryonic 
Newtonian acceleration scales within the DM scale radius by considering a large sample 
of late -type galaxies. We rely on a Markov Chain Monte Carlo (MCMC) method to 
estimate the parameters of the halo model and the stellar mass - to - light ratio and 
then propagate the uncertainties from the rotation curve data to the estimate of the 
acceleration scales. This procedure allows us to compile a catalog of 58 objects with 
estimated values of the B band absolute magnitude Mb, the virial mass M V i r , the DM 
and baryonic Newtonian accelerations (denoted as gDAi(ro) and gbar(ro), respectively) 
within the scale radius rn which we use to investigate whether it is possible to define 
a universal acceleration scale. We find a weak but statistically meaningful correlation 
with M V i r thus making us argue against the universality of the acceleration scales. 
However, the results somewhat depend on the sample adopted so that a careful analysis 
of selection effects should be carried out before any definitive conclusion can be drawn. 
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1 INTRODUCTION 

Galactic systems cover a wide range in size, mass, lumi- 
nosity, colours and morphology so that it is almost obvi- 
1 ous that finding a single model able to fully reproduce this 
, rich phenomenology is doubtful. It is nevertheless some- 
what surprising that some order in this wide parameter 
space may be created thanks to the existence of statisti- 
cally meaningful correlations among galaxy structural and 
dynamical parameters. Textbook examples are represented 
by the Faber - Jackson (1976) and Tully- Fisher (1977) re- 
lations between velocity indicator and luminosity for early 
and late - type galaxie s, respectively, the Fundame ntal Plane 
|Dressler et al. I Il987i ; iDiorgovski fc Davis 1 1 19871 ) between 
velocity dispers ion, size and effective intensity and the pho- 
tometric plane (|Khosroshahi et al. l2000l ; lGraham l2002t ) be- 
tween photometric parameters of elliptical systems. Besides 
their usefulness as possible distance indicators, the existence 
of scaling relations among galaxy structural parameters is 
evidence that some hitherto not fully explained mechanism 
has taken place during the formation and evolutionary pro- 
cesses leading to the currently observed correlations. Find- 
ing and interpreting scaling relations has therefore recently 
become an active research field in the literature of galaxy 
properties. 

* Corresponding author : winnyenodrac@gmail . com 



Motivated by these considerations, 

iKormendv fc Freeman I l|2004l ) described the dark halo den- 
sity profile pDM of Sc - Im and dwarf spheroidal galaxies with 
the pseudo- isothermal profile, pdm{t) = pQr%j(r 2 + r 2 .), 
to fit the rotation curve data of 55 systems. Under the 
maximal disk hypothesis, they found that the quantity 
po = poVo, proportional to the halo central density, stays 
remarkably constant at the value po — 100 Mq/pc 2 for all 
the galaxies in their sample. Donato et al. (2009, hereafter 
D09) then further extended this result considering the 
coadded rotation curves of ~ 1000 spiral galaxies, indi- 
vidual mass models for spiral and dwarf systems and the 
w eak lensin g data of a set of elliptical galaxies. Adopting 
a iBurkert I (| 19951 ) model for the dark halo, they found 
Po r o = 141^52 Mq/pc 2 over 14 units of B magnitude (i.e., 
14/2.5 = 5.6 orders of magnitude in luminosity) where it 
is important to stress that (po,ro) are defined with respect 
to the Burkert rather than pseudo- isothermal profile. The 
D09 result may also be expressed in terms of the Newtonian 
acceleration at the scale radius ro . This is given by : 

Jm/ r = s-i — - = 3.21 j 2 x 10 cm/s 

and should stand out as a universal quantity provided that 
D09 result is correct. This result has since been generalized 
to the gravitational acceleration of the baryons within ro 
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by Gentile et al. (2009, hereafter G09). Using a sample of 
spiral galaxies with carefully measured rotation curves and 
excluding some systems based on quality criteria, G09 found 



g ba r(r ) = 5.7li;i x 10 



-10 / 2 

cm/s 



to be constant for all the galaxies in their sample. 

As already stated above, both the D09 and G09 results 
were obtained adopting the Burkert density profile, 



P{r) = 



Porl 



(r + r )(r 2 + r 2 ) 



for the DM haloes. According to D09, this choice is mo- 
tivated by the conducive property of the Burkert profile 
to mimic a pseudo - isothermal core in the inner regions 
and a r -3 scaling in the outer ones, s imilar to the popu- 
lar NFW profile (jNavarro et al. I ll997l V Moreover, such a 
model is k nown to make extremely good fits to rota tion 
curve data (jGentile et al. Il2004 [20071 ; iDel Popolo 1120091) al- 
though some few remar kable exceptions exist ( Simon et al. I 
120051 : IDel Popolo II2012I ) . 

The D09 and G09 results have also been revised and 
criticized by different authors. First. iBovarskv et al. I (|2009l ) 
extended the analysis to group and cluster scale systems 
evaluating the column density defined as : 



sible their original analysis, but improving some substantial 
points. First, we almost double the G09 sample which also 
allows us to investigate possible selection effects. Second, we 
estimate the Burkert halo model parameters by directly fit- 
ting the rotation curve data rather than adjusting previous 
results already available in the literature and based on the 
use of different density profiles. Third, we adopt a Bayesian 
fitting procedure to infer a realistic estimate of the errors on 
the quantities of interest also taking into account the uncer- 
tainties on the mass - to - light ratio. Rather than considering 
a column density, we express our results in terms of the New- 
tonian acceleration at the DM scale radius and investigate 
whether this quantity correlates with the luminosity, mass 
and size of the galaxies. 

The plan of the paper is as follows. Sect. 2 gives some 
general consideration on how we define and compute the 
Newtonian acceleration also explaining to what extent this 
quantity may or may not be considered model independent. 
We then present the rotation curve sample in Sect. 3 and 
discuss the results of fitting our adopted galaxy model (in- 
cluding gas, disc and DM components) to the data. Scaling 
relations are investigated in Sect. 4, while our summary and 
discussion are finally presented in Sect. 5. In Appendix A, 
we list the constraints on the model parameters and the es- 
timate of the Newtonian accelerations so that the interested 
reader can use them for any further investigation. 
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which is proportional to the mean surface density within r, 
(given by p*r*) and exactly equal to po?"o for the Burkert 
model. A fit to their sample data gives 

with S in M©/pc 2 thus contr adicting the D09 results o n the 
constancy of po r o- Similarly. iNapolitano et al. I (|20 lOf ) con- 
sidered a large sample of elliptical galaxies and estimated 
the quantity (ppM )R e f f which is proporti onal to the column 
density defined bv lBovarskv et all |2009l '). but using the ef- 
fective radius Reff as reference scale. Although they agree 
with D09 that this quantity is constant for spirals, these au- 
thors nevertheless found that early - type galaxies violate the 
constant density scenario by a factor ~ 10 on average and 
by a factor ~ 5 if one only considers systems spanning the 
sa me mass regime. A further analysis has been carried out 
bv lCardone fc Tortoral I^OIOT ) based on strong lensing and 
velocity dispersion data for intermediate redshift lens galax- 
ies. They found that the column density correlates with the 
mass and the luminosity although the slope of the scaling 
relations depend on both the DM halo model adopted and 
the stellar IMF. 

The original D09 and G09 results and the following tests 
quoted above can not be straightforwardly compared be- 
cause of the differences in the reference radii adopted, the 
halo model and the mass and luminosity ranges covered. As 
such, drawing a definitive conclusion on the existence or not 
of a universal quantity is actually not an easy task. In an at- 
tempt to give a clear cut answer to this question, we revisit 
here the D09 and G09 results by mimicking as close as pos- 



2 NEWTONIAN ACCELERATION SCALES 

The concept of Newtonian acceleration is so well known that 
spending more than few words on it could seem quite su- 
perfluous. Nevertheless, while its definition is rather sim- 
ple, some subtleties arise when dealing with actual galaxies. 
First, spiral galaxies may be considered three component 
systems with stars and gas representing the baryon mass 
and the halo accounting for the dark matter term. In a first 
and reasonably good approximation, stars are distributed 
in a thin disc well modeled by the Freeman (1970) exponen- 
tial law, while the gas component roughly follows the disc 
but with strong irregularities eventually leading to a clumpy 
distribution. Modeling such an irregular profile with an an- 
alytical expression could introduce a bias in the derivation 
of the Newtonian acceleration thus altering the analysis of 
the scaling relations we are interested in. We therefore use a 
different approach noting that, for a spherically symmetric 
system, the circular velocity simply reads 

vl{r) = GM(r)/r . 

Since observations of the luminosity profile directly provide 
us the circular velocity of both baryonic components, we can 
define an effective mass as 

M eff (r)=rv 2 c (r)/G 
and then estimate the corresponding acceleration as 

g(r) = GM eff (r)/r 2 . 
We therefore define the Newtonian accelerations as : 



DM and stellar accelerations 3 




Figure 1. Dark matter acceleration ratio for the Burkcrt vs NFW models assuming the same virial velocity. Short dashed, solid and 
long dashed lines refer to models with c/vfw = 2.5, 5.0, 7.5, respectively, while we set C2{Bur) / c 2 {N FW) = 0.5, 1.0, 2.0 from left to right 
panels. Dotted lines delineate the region where the ratio is equal to unity within a 10% confidence. 
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M DM (r) = M vir x 



(2) 



with M V i r and R v i r = [(3M v i r )/(4nA V i r p cr it)] 1 ^ 3 the virial 
mass and radiufl respectively, and 



where Vdisc and v gaB are the measured stellar disc and gas 
circular velocities at the radial distance r from the galaxy 
centre. It is worth stressing that Eqs.fTJ) are exact only for 
a spherical mass distribution, while neither the disc nor the 
gas are spherical systems. Although not being formally cor- 
rect, such definitions possess the remarkable property to be 
directly related to the data thus avoiding the errors due 
to an incorrect modeling of the mass distribution related 
to gradients in the stellar M/L ratio or dumpiness in the 
gas profile. We nevertheless stress that, when fitting rotation 
curve data, we do not rely on the approximated Eqs.([T]), 
but rather correctly compute the theoretical circular veloc- 
ity from the gravitational potential obtained by solving the 
Poisson equation for the given componenlQ. Should we have 
used the spherical approximation, significant biases could 
have been introduced in the estimate of the halo parameters 
depending on the details of the disc and gas distribution 
properties (such as, e.g., the disc scalelength and the gas 
dumpiness). 

Computing the Newtonian acceleration for the dark 
matter term presents a different problem. On the one hand, 
it is customary to assume that the halo is spherically 
symmetric so that the acceleration is simply qdm(t) = 
GMdm(j) /r 2 with Mdm(?) the halo mass within r. On the 
other hand, the halo circular velocity is not directly mea- 
sured so that one can not resort to the analog of Eq.{TJ to 
infer the DM acceleration from the data. One is therefore 
forced to assume a halo model, determine its parameters by 
fitting to the observed total rotation curve and then esti- 
mate the accelera tion. In o r der to be consistent with G09, 
we will assume a iBurkert I l|l995l ) model so that the mass 
profile reads : 



H(r) = In H 

ro 



r \ 1 
-arctan | — +- In 
.r J 2 



1+ ( - 

ro 
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with ro a characteristic radius. The model is assigned by 
the two parameters (M v i r , ro), but, when fitting to the rota- 
tion curve data, it is more convenient to reparameterize the 
model in terms of the circular velocity at the virial radius, 
V 2 ir = GM vir / R V ir, and the concentration c vir = R vir /ro- 
The DM Newtonian acceleration can then be written as 



gDM(r) = 



GM vir fi(r) 



Qvir X 



H(r)/r 



(4) 



r 2 ^(R vlr ) au " " ^R m r)/R 2 vir 

having set g v i r = GM v i r / ' R 2 ir = V 2 ir / 'R V ir for the acceler- 
ation at the virial radius. Eq.@ allows us to immediately 
estimate how the acceleration depends on the halo model. 
First, we note that, since the virial velocity V v %r is typically 
well determined by fitting the rotation curve, its value can 
be taken as model independent. That is to say, in order to 
have two different models fitting the same data, they must 
have the same virial velocity. We will therefore consider two 
different models with the same V V i r (and hence the same 
R V i r and g V i r values), but two different mass profiles de- 
scribed by the functions /ii(r) and ^(r). The ratio among 
the corresponding accelerations will then be : 



flolfM = Hl(r)/Hl(Rvir) 
SdmM ^2{r) / ^ 2 {Rvir) 



(5) 



As an example, we consider here the Burkert and NFW 
jNavarro et al. |[l997l ) as models 1 and 2, respectively, thus 



setting 



Actually, for most cases, the original references for the data 
already provide also the circular velocity of the disc and gas com- 
ponents assuming a fiducial stellar M/L ratio T so that we have 
only to rescale them to match our estimated T. 



/te(r) = Inll + — --p [ 1 + —r ( 6 ) 

2 We follow [Bryan &; Norman to set A^jy, while p C rit — 

3Hq/8ttG is the critical density of the universe. 
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with cnfw = Rvir/r s and r s the characteristic radius of the 
NFW profile. In order to meaningfully compare the models, 
we first introduce a new concentration parameter defined 
as C2 = R v ir/R2 with Ri the radius where the logarithmic 
slope of the density profile, a = — dlnp/dlnr, equals the 
isothermal value a — —2. It is easy to check that 



Rvir 

~R^ 



= < 



Rvir Cvir 



1.52r 1.52 

Rvi 



for Burkert 



cnfw for NFW 



Fig. \T\ shows the ratio S^mM /SamM f° r the Burkert vs 
NFW models as a function of r/R v i r assuming the two 
models have the same virial velocity and different values 
of the NFW concentration cnfw and C2{Bur) /c2(NFW). 
This plot highlights an important feature of the Newtonian 
acceleration, namely its dependence on the particular halo 
model adopted. This point clearly has to be taken into ac- 
count when looking for scaling relations among Qdm{t) and 
stellar or DM quantities. The central panel, however, sug- 
gests a possible way out of this problem. Indeed, if we ask 
that two different models fit the same rotation curve data, 
their circular velocity profiles must be similar, at least over 
the radial range probed by the data. We have checked by 
trial and error that, for the Burkert and NFW models, this 
is the same as asking that the virial velocity Voir is the same 
and the C2 concentration is comparable. In such a case, the 
central panel in Fig.[T]shows that gfjM ( r ) — QdaF^ ( r ) within 
10% if r > r min with r min ~ 0.1 — Q.2R v i r (depending on 
the cnfw value). In the following, we will use the Burkert 
model to estimate the acceleration values, but will investi- 
gate scaling relations only for values of r > r m in in order 
to be confident that a mismatch in the modeling does not 
significantly bias the results. 



3 ESTIMATING THE ACCELERATION 

Both the baryon and DM accelerations require rotation 
curve data in order to be evaluated. Indeed, HI observations 
directly give the amount of gas and its profile thus allowing 
us to infer v gas (r), while photometry and an estimate of the 
stellar M/L ratio gives Vdisc(r). Adding in quadrature these 
terms to the theoretically predicted DM circular velocity al- 
lows us to get the combined circular velocity to be compared 
with the rotation curve and hence determine the halo model 
parameters (Voir, Cmr)- In the following, we will first describe 
the galaxy sample we used and the fit procedure and then 
will explain how we estimated the acceleration scales. 



bias the estimate of the halo parameter^ . We finally end up 
with 58 galaxies from different sources (Ide Blok fc Bosma 



20021 ; [Simon et al. 1120051; ISpano et al I|200S| ; Ide Blok et al 



20081 ; lOh et al. Il201ll ; ISwaters et al. M2OIII ) spanning nine 
units of B band magnitude (roughly, —22 < Mb ^ — 13) 
and 9/2.5 orders in luminosity and different spiral classes 
(from dwarfs system to low and high surface brightness 
galaxies) . 

It is worth stressing that our compilation is not based 
on any selection criteria other than the availability of mea- 
surements covering a large radial range. One could, however, 
wonder whether such a blind selection may lead to reliable 
results. Indeed, poor quality curves may be fitted by models 
with unreliable parameters so that the estimate of the ac- 
celeration scales should not be trusted. On the other hand, 
selecting a priori only a subset of the galaxies may risk bias 
in the analysis of scaling relations excluding systems which 
drive an eventual weak correlation. In order to investigate 
this issue, we will start by retaining all the galaxies and then 
look for the impact of different quality cuts on the slope and 
scatter of the scaling relations among acceleration scales and 
other parameters of interest. 

In order to determine the model parameters, we 
carry out a likelihood analysis by maximizing C(p) oc 
cxp[- X 2 (p)/2] with 



AT 



s (r«) 



(7) 



where v° bs (ri) is the measured circular velocity for the i-th 
point (with Oi the error), vl h (ri,p) the theoretically pre- 
dicted value for the given set of parameters p and the sum 
is over the Af data points. It is worth stressing that, since the 
errors reported from different authors have been obtained in 
different ways and typically take care of possible systematic 
effects, the reduced \ 2 f° r the best fit model is not expected 
to be unity. 

The fitting procedure allows us to determine the halo 
model parameters (V v i r , ew) and the stellar M/L ratio T*. 
This could seem at odds with our previous statement that 
the disc circular velocity is directly estimated from the data. 
Actually, this is only approximately true. Indeed, the pho- 
tometry of the galaxies allows to set the disc profile and 
hence evaluate the corresponding circular velocity provided 
T* is known. This latter quantity can still be guessed from 
the multiband photometric data using a prescription to con- 
vert colors into a fiducial T*. However, uncertainties in the 
properties of the stellar population (such as, e.g., its age 
and metallicity) makes the estimated T* affected by a large 
error. To take care of this issue, rotation curve papers typ- 
ically report the values of Vdisc(r) for T* = 1 so that the 
actual stellar M/L ratio has to be determined by the fit to 



3.1 Data and fitting method 

We have searched the literature for systems with high qual- 
ity rotation curve data probing the gravitational poten- 
tial with good sampling, large radial extent (i.e., up to 
R > Ropt = 3.2Rd with Rd the disk scale length) and small 
errors. Whenever possible, we rely on Ha data or a combi- 
nation of HI and Ha so that the potential impact of beam 
smearing on the circular velocity is reduced and does not 



3 A possible concern on the use of the Ha data is related to Ha 
possibly tracing biased star forming regions of the disc (as, e.g., 
spiral arms). As a consequence, the measured circular velocity 
could be non representative of the disc stars kinematics. While 
we can not definitively exclude such a possibility, it is nevertheless 
worth noting that, whenever both HI and Ho data are available 
in the same region, they closely track each other. We are therefore 
confident that no significant bias is induced by the use of Ha data 



DM and stellar accelerations 5 



NGC3621 NGC5055 NGC 925 




5 10 15 20 25 10 20 30 40 2 4 6 8 10 

R (kpc) R (kpc) R (kpc) 



Figure 2. Best fit models superimposed on the data for three example galaxies. Red solid, blue long dashed, purple dot dashed, black 
dotted and green short dashed lines refer to the the full, disc, gas, halo and bulge contributions to the the circular velocity. 



the data. We therefore end up with a three parameters set, 
i.e., in Eq.j7|), it is p = (Y*, Kir, <w). 

In order to efficiently explore this 3D space, we use 
a Markov Chain Monte Carlo (MCMC) method run- 
ning three chains check ing the convergence through the 
iGelman fc Rubin! l) 19921) criterion. After cutting the chains 
to avoid the burn - in phase and thinning them to erase spu- 
rious correlations, we merge the three chains in a single one 
which can be used to infer constraints on the model pa- 
rameters. Note that, in a Bayesian framework, the best fit 
parameters do not represent the most reliable estimate 
of the single parameters pi. On the contrary, one should look 
at the marginalized likelihoods obtained by integrating £(p) 
over all the parameters but the i - th one. Using the MCMC 
method, this reduces to look at the histogram of the i-th 
parameter values to infer the median and the 68 and 95% 
confidence ranges. Although we will report also the best fit, 
in the following, we will consider the median and 68% confi- 
dence range as our final estimates of the model parameters. 

3.2 Results from rotation curve fitting 

It has long been known that cored models are quite efficient 
at fitting spiral galaxies rotation curves. It is therefore not 
surprising that the adopted Burkert profile provides a good 
fit for almost all of the galaxies in our sample. It is neverthe- 
less worth spending some more words on how one can decide 
whether the match between the data and the best fit model 
can be considered as good or not. As already stressed above, 
the errors on the measured circular velocity are not Gaus- 
sian distributed since they can include contributions for the 
asymmetry between the approaching and receding sides of 
the galaxy. As a first consequence, one can not rely on the 
usual reduced x 2 criterion (i.e., \ 2 — X 2 /d.o.f. ~ 1) to deem 
a fit as good or bad. Moreover, our sample is comprised of 
galaxies collected from different authors, each one of them 
using their own recipe for assigning errors to the rotation 
curve data. As such, it is not straightforward to work out a 
single quantitative criterion which works equally well for all 
cases. We have therefore relied on three different qualitative 
criteria to deem the matching among data and best fit model 
as evidence in favour of the Burkert profile. First, we look 
at the x 2 value. Indeed, although its normalization depends 
on how the errors have been evaluated, it is nevertheless 
clear that a large value is indicative of a poor fit (unless the 
errors have been grossly underestimated which is not the 



case). A small x 2 value can nevertheless be obtained if the 
errors have been overestimated. As a consistency check, we 
will therefore look at the root mean square of the residu- 
als and their trend with the radius. For a model well fitting 
the data, we indeed expect that the residuals (scaled with 
respect to the measurement errors or to the observed circu- 
lar velocity values) are small and scatter around the mean 
value. As such, the rms of the residuals will be small and no 
correlation with r can be observed. On the contrary, a trend 
with r will be obtained if the model is systematically failing 
to fit the data. 

A few examples will help to clarify the use of the above 
quality criteria. The left panel in Fig. [2] shows the rota- 
tion curve data for NGC 3621 and the superimposed best 
fit model (x 2 = 0.75). This is a typical case where one could 
rely on the expectation based on the reduced \ 2 value are 
in accordance with our definition of well fitted system. The 
best fit model closely tracks the data and there is no trend 
of the very small residuals with the distance from the galaxy 
centre. On the contrary, the case of NGC 5055 (shown in the 
central panel of Fig. [2]) is an example of a system with a large 
reduced \ 2 (x. 2 = 3.63), but with a nevertheless reasonably 
well fitted rotation curve. Indeed, the best fit model tracks 
the data almost everywhere but in the very inner regions 
where the circular velocity is likely to be mainly determined 
by the bulge component which we hold fixed in our analy- 
sis. We have therefore included this galaxy in our sample of 
well fitted systems since the divergences may be very likely 
an evidence of a failure in the bulge rather than the halo 
modeling. Finally, the right panel of Fig. [2] shows the data 
and the best fit model for NGC 925 which is a good example 
of why one must not only rely on x 2 value. Indeed, we get 
X 2 = 1.31 for this galaxy (smaller than for NGC 5055), but 
the best fit model fails to follow the data overestimating the 
circular velocity in the inner regions and underestimating 
it in the outer ones. Moreover, the best fit circular velocity 
profile keeps almost linearly increasing, whereas the data 
shows some hints of a flattening of the rotation curve. Mo- 
tivated by these discrepancies, we therefore consider this fit 
an example case of a poorly fitted system. 

Applying these conservative criteria, we finally find that 
all the galaxies but six are well fitted by the Burkert profile. 
This is in agreement with the previous results in the litera- 
ture which have shown that this halo model does a very good 
job in fitting the rotation curves of spiral galaxies of different 
morphological classes. Although this could be considered as 
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further evidence in favour of cored vs cuspy models, it is 
worth stressing that tackling this unresolved debate is not 
our aim here. From our point of view, having successfully fit- 
ted the rotation curve data of such a large sample of spiral 
galaxies allows us to confidently rely on the Burkert pro- 
file to estimate the Newtonian acceleration scales of interest 
here. 

3.3 Acceleration scales 

Having constrained the disc stellar M/L ratio and the 
halo model parameters, we can now use Eqs.JT} and Q 
to estimate the Newtonian acceleration of both the baryon 
(disc + gas) and DM components. To fully take into account 
the uncertainties on the fitting parameters and the corre- 
lations among them, we evaluate gi(rj) along the Markov 
chain for each galaxy and use the resulting values to es- 
timate the median and 68% confidence ranges. Note that, 
here, the label i denotes the galaxy component (baryons or 
DM), while j refers to the reference radius adopted. Two 
choices are possible. First, one can set rj = Rd, i.e. esti- 
mate the acceleration at the disc scalelength radius. On one 
hand, such a choice could be well motivated noting that 
the data mainly sample the intermediate disc regions (i.e., 
R ~ few Rd) so that gi(Rd) is likely to be well constrained. 
Actually, this is also the region where uncertainties on the 
baryon components (such as the disc stellar M/L and the 
bulge modeling) have the larger impact. We therefore pre- 
fer to set r = ro in order to be consistent with the G09 
analysis. It is worth stressing, however, that, for the typi- 
cal concentration values we have obtained from the fitting 
analysis, one gets ro/R V i r < 0.1 so that one is actually esti- 
mating the acceleration in the region which is most sensitive 
to changes in the halo model as can be seen from Fig.Q] We 
therefore warn the reader to not rely on any scaling rule to 
infer the gDM(ro) values for a different model from the ones 
we find here using the Burkert profile. 

Although the data on the rotation curve typically ex- 
tend up to large r/Rd values, it is nevertheless possible that 
the constrained ro is larger than R ma x, this latter quan- 
tity being the last measured point of the rotation curve. 
In such a cases, one could still rely on the estimated ro to 
infer the <?i(ro) values, but this comes at the price of assum- 
ing that the model fitted up to R max can be extrapolated 
outside. In order to avoid any possible systematic error re- 
lated to such a choice, we prefer to exclude the galaxies with 
ro > Rmax even if their rotation curve are well fitted by the 
assumed model. For completeness, we give in Table, [A"2l the 
constraints on (/i(ro) for all the galaxies in our sample, but 
we will use in the rest of the paper only those systems hav- 
ing both a well fitted rotation curve and ro < Rmax referred 
to in the following as the Good sample. 



gDM(ro) = 3.2l}'2 x 10 -9 cm/s 2 as inferred from a sam- 
ple of 28 spiral galaxies along the Hubble sequence. Using 
our larger sample, we can revisit their results, but rely on a 
somewhat different procedure. First, we note that G09 did 
actually not perform any fit to the rotation curve data, but 
relied on the literature and an approximate law to convert 
the scale length of the model used by other authors to their 
Burkert model ro radius. While such a scaling law gives a 
correct local matching among different models, it is worth 
stressing that this does not automatically ensure that the 
theoretically predicted rotation curves well match also in 
the outer regions. On the contrary, here, we fit the full ro- 
tation curve with the Burkert model thus being also able 
to check whether or not the model well reproduces the ob- 
served data. As a result, we find that the scaling law used by 
G09 indeed provides values of gi(ro) consistent with our ones 
with (A log 3dm (ro)) = -0.03 and (Alogflw(r )) = -0.02 
as estimated by the 18 galaxies common to both samplefl 
However, in order for such an agreement to hold, one has 
first to check that the fit to the data is good, which is not 
the case for the full sample. As such, some of the galaxies 
(see TabIe [A"2|) included by G09 must actually be excluded 
in order not to bias the final results. 

A further difference with G09 concerns the estimate of 
the uncertainties. Since we have determined rather than as- 
sumed the stellar M/L ratio and the halo model parameters, 
we have been able to attach to each (/i(ro) value an error fully 
taking into account both the fitting parameters uncertain- 
ties and their correlations. In order to be as conservative as 
possible, we then add in quadrature a systematic error re- 
lated to the uncertainties on the method used to determine 
<7i(ro). Such a systematic is hard to ascertain so that we 
empirically rely on the comparison with the G09 values for 
the galaxies common to both samples and add in quadrature 
0.23 to the statistical errors on log <?i(ro) since this is the mis 
of both components of Alogpi(ro). The final uncertainties 
turns out to be smaller than the ones adopted by G09 based 
on a qualitative analysis of the impact of the M/L and sys- 
tematic uncertainties. We are, however, confident that our 
estimate of the errors is better motivated by a statistical 
point of view so that we prefer to rely on our guess in or- 
der not to underate the ability of the data to constrain the 
variation of gi(ro). 

Should the acceleration scales gi(ro) be universal quan- 
tities as argued by G09, one should find no correlation with 
any stellar or DM property. In order to check their hypoth- 
esis, we therefore fit the following log - linear relations 

log 5i (ro) = alog( I5 ^|-)+6 1 
log 3i(r ) = aM B + b , 



4 SCALING RELATIONS 

The analysis first laid out in G09 (but see also 
iKormendv fc Freemanl |2004| ; iDonato et all 120091 for pre- 
liminary investigations) has concluded that both the baryon 
and the DM accelerations evaluated at the Burkert scale 
radius ro stay constant over almost 14 B magnitudes in 
luminosity with gbar(ro) = x 10 -10 cm/s 2 and 



4 We overlook 10 galaxies from the G09 sample because of dif- 
ferent reasons. First, we do not use the six Milky Way dwarf 
satellites since the estimate of <7i(ro) has been obtained relying 
on a fit to the velocity dispersion data, while here we are fitting 
rotation curves which are not available for these systems. We also 
reject four more galaxies since we have been unable to retrieve 
the rotation curve data so that we can not estimate the stellar 
M/L ratio and the halo model parameters. 
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Figure 3. Best fit models superimposed on the data for the DM acceleration gDAt( r o) as a function of the halo virial mass logM„; r n = 
log (M v i r /10 Mq) and the galaxy absolute B band magnitude Mj. The solid (dashed) line refers to the Baycsian (direct) fit. 




l0 8 M .ir.ll M„ ~Z 

Figure 4. Left and centre. Same as Fig. [3] but for the baryon acceleration gbar{ r o)- Right. Our best fit gt, ar (ro) vs Mb relation 
superimposed to the G09 (red) and our (black) acceleration data. 



with i — {DM, bar} using the Bayesian methocjfl in 
ID' Agostini I HH). Note that such an approach allows us 
to take care of the uncertainties on both fitting variables 
and also to estimate the intrinsic scatter of the correlation 
in an unbiased way. The results are summarized in Table[T] 
for three different galaxy subsamplcs. 

Let us first consider the Good sample made out of all 
the galaxies well fitted by the Burkert profile and with a 
scale radius ro smaller than the radius of the last measured 
point on the rotation curve. As it is apparent from both the 
slope values and Fig.[j|] the hypothesis that the DM accel- 
eration gDM(ro) is a universal quantity is not supported by 
the data. Indeed, from a statistical point of view, the ab- 
sence of any correlation (whether be it with the halo mass 
M V ir or the galaxy luminosity Mb) can be safely excluded at 
the 95% CL. Such a result could seem in fierce disagreement 
with the G09 one. Actually, this is not the case. Indeed, 
fitting their smaller galaxy sample, G09 obtained a slope 
a ~ 0.0 ± 0.2 for the <?dm(^o) vs Mb so that, although the 
median value is different, the two results are statistically 
consistent. Moreover, the intercept of our gDM(ro) vs Mb 



5 As a consistency check, we also use the direct fit method, i.e., we 
minimize the usual \ 2 neglecting the uncertainties on the fitted 
quantities and the intrinsic scatter. Since the results on the fit 
parameters (a, b) turn out to be in good agreement, we will only 
refer in the following to the results from the Bayesian method. 
We nevertheless show both best fit models in the relevant figures. 



correlation translates into <7dm(»"o) = 5.OI24 x 10 -10 cm/s 2 
which is well consistent with the claimed constant value 
found by G09 and quoted before. Given the almost flat slope 
and the smaller sample used by G09, it is not surprising that 
their median value for a shifts towards a null one so that 
we argue that there is no significant difference among the 
two results of the (?da/(Vo) vs Mb relation. What is, nev- 
ertheless, different is their interpretation. Indeed, the weak 
correlation between (?da/(Vo) an d Mb can be qualitatively 
considered as a lack of correlation thus motivating the status 
of a universal quantity given to <?dm(»"o) in G09. However, 
the significant correlation with the halo virial mass can be 
taken as strong evidence against the universality of gDM(ro) 
even if one neglects the very weak but yet statistically mean- 
ingful scaling with Mb- We therefore conclude that no DM 
universal acceleration scale may be defined. 

While we substantially agree with G09 concerning the 
scaling of gDAi(ro) with Mb, there is on the contrary a 
rather strong disagreement for the scaling of the baryon 
acceleration gtar{ro). In their paper, G09 claim that this 
quantity is constant over 14 units of B mag, while Fig.0] 
and the slopes in Table[T] make us definitely argue against 
this conclusion. Indeed, a vanishing slope is well outside the 
95% CL for both the gb ar (ro) - M vir and gb ar {ro) - Mb rela- 
tions. Moreover, the a values are not close to zero so that 
taking a constant gw(ro) leads to a severe mismatch to the 
data. 

It is worth wondering why we find a result inconsis- 
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Table 1. Scaling law parameters as output from the fit to three different galaxy subsamples. For each fit, we give the median and 68% 
confidence range for the slope a, the intercept b and the intrinsic scatter er^t obtained using the D' Agostini (2005) Bayesian method. 
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Figure 5. A fundamental plane -like relation for <?i(ro) as function of Mb and R^. 



tent with the G09 ones. To this end, we reduce our sample 
by only including the 16 galaxies common to both catalogs 
and fit the scaling relations to this subset. As a result, we 
find that the slopes of the <?dm(^o) - M„i r and Qdm{to) - Mb 
relations are now significantly smaller and consistent with 
null values within the 95% CL so that we recover the G09 
result that this quantity can be considered as roughly uni- 
versal. However, although the scaling laws become shallower, 
gbar{ro) still significantly correlates with both the halo virial 
mass and the galaxy luminosity so that the assumption of 
constant baryon acceleration can not be retained. 

In order to trace the origin of this difference, it is worth 
first remembering that the G09 sample has been assembled 
starting from a larger catalog and cutting out objects which 
do not fulfil some quality selection criteria. On the contrary, 
the only criteria adopted to include a galaxy in our Good 
sample relies on the goodness of the fit with the Burkert 
profile. Should the selection of galaxies be responsible for 
the different slopes, the fit to the 16 galaxies in common 
with the two samples would have returned consistent slopes 
which is not always the case. To further explore this issue, 
we have plotted in the right panel of Fig.[4]the G09 gbar(ro) 
values superimposed on our dat£0. An interesting lesson can 
be drawn from this figure. Should we exclude the galaxies 



6 This plot also shows us whether the estimated gi,ar( r o) values 
for the galaxies in common to both the G09 and our Good sample 
agree or not. In the former case, two almost superimposed point 
will appear in the plot. 



with log gbar(ro) > —8.0 and add the Milky Way dwarfs with 
Mb > —12, a constant gbar{ro) would be preferred. All the 
four systems with log gbar{ro) > —8.0 (namely, NGC2841, 
NGC3521, NGC4736, NGC6946) belong to the THINGS 
sample so that they have a high quality rotation curve and 
a not small inclination (as estimated from the HI data). It is, 
however, nevertheless possible that the presence of warps or 
non circular motions due to, e.g., bars and/or other sources 
of instabilities make the rotation curve (notwithstanding the 
quality of the data) a poor estimator of the underlying gravi- 
tational potential thus motivating the exclusion of these sys- 
tems from the G09 sample. In the range —20 ^ Mb ^ —15, 
our data are well superimposed on the G09 ones except for 
few cases where our estimates are larger than those of G09. 
However, since the values in G09 comes from scaling the re- 
sult for other halo models to the Burkert one, while our own 
are directly based on fitting this halo profile to the rota- 
tion curve data, we consider our values more reliable. Based 
on such a comparison, we argue that the most likely reason 
for the discrepancy on the constancy of gbar{ro) with Mb 
is mainly due to our sample including galaxies with larger 
acceleration values. However, in order to definitely discrimi- 
nate between the two contrasting results, one should comple- 
ment our sample with systems probing the luminosity range 
Mb > —12 to test whether the inferred gbar(ro) — Mb scal- 
ing can be extrapolated to this regime or a constant value 
for the baryon acceleration scale should be preferred. It is 
also worth remembering that B band is a poor estimator of 
the stellar mass M* which is actually a far better indicator 
of the galaxy status. Being M* hard to estimate in a model 
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independent way, one should better rely whenever possible 
on Spitzer 3.6/im. magnitudes which are a better tracer of 
the stellar mass. 

As a further test, we have investigated whether the in- 
clusion of systems with poorly determined accelerations can 
bias the inferred constraints on the scaling relations. Actu- 
ally, the S/N ratio is typically quite large so that there are no 
systems prone to large errors. We have nevertheless selected 
a subsample (referred to as our Gold sample) requiring the 
S/N ratio to be larger than 50 for both gi(ro). Repeating the 
different fits for this subsample, we find the results in Ta- 
ble]]] which are in good agreement with those from the Good 
sample. We will not discuss further this case in the follow- 
ing, but here we specifically note that the intrinsic scatter 
decreases for all the correlations investigated. 

The scaling relations we have found above are affected 
by a significant scatter. In order to see whether it is possible 
to reduce it, we look for a 3D correlation among the accel- 
eration scales and the two stellar parameters directly mea- 
surable from the galaxy photometry, namely the total lumi- 
nosity Mb and the disc scale length radius . We therefore 
use the same Bayesian fitting method to determine the co- 
efficients of the log - linear relations 

log ~. (r ) = aM B + b log R d + c . 
Using the Good sample, we find (median and 68% range) : 



-o.o88±g:gS , b - 



c= -9.794+°:™ 



for i = DM and 



n 071 +0.038 
— U.Z, I 1_().040 



c = -i3.4Big:2S 



fl D79+ - 038 
-U.UfZ_Q.Q42 i 



— U.114_ .058 



n icq+0.066 
~U.iO__o.062 



= 0.246] 



for i = bar (see Fig. [5]). The addition of a second parame- 
ter reduces the intrinsic scatter for both <?DA~(ro) and (to 
a lesser extent) gbar{ro), but one has also to take into ac- 
count that a smaller number of galaxies has been used since 
Rd is not available for all the systems in the Good sample. 
A definitive conclusion on whether a 2D or 3D scaling law 
is better suited to forecast the acceleration scales could be 
possible only increasing the statistics and carefully modeling 
the uncertainties on the quantities involved. 



5 CONCLUSIONS 

Scaling relations among dark matter and stellar quantities 
can offer valuable insights into the galaxy formation process. 
Similarly, the existence of universal quantities may be taken 
as evidence of some hitherto unknown underlying physical 
mechanism coupling the evolution of the dark matter and 
baryons in the galaxy. Motivated by the results of D09 and 
G09 and given the presently confused situation, we have here 
revisited the issue of the universality of DM and baryon 
Newtonian acceleration evaluated at the DM scale radius. 



For this purpose, we have used a large sample of late - type 
galaxies spanning roughly nine orders of B magnitude. The 
use of a MCMC procedure has allowed us to both estimate 
the model parameters (namely, the stellar M/L ratio and the 
DM halo virial mass and concentration) and then correctly 
propagate the errors from both the rotation curve and model 
parameters to the final estimates of the acceleration scales 
gi( r o) (with i = DM, bar for DM and baryons, respectively). 

Using a sample made out of well fitted rotation curve 
galaxies, we find that gDM{ro) weakly correlates with the B 
band absolute magnitude Mb in satisfactorily good agree- 
ment with what is expected from the claimed constancy of 
pofo found by D09. However, we also find a strong correla- 
tion with the virial mass M V i r which makes us argue against 
the universality of _£.__(»_)■ The slopes of the gbar(ro) vs 
Mb and M V i r relations turn out to be definitely non van- 
ishing so that we can exclude at 95% CL that this quantity 
is a universal constant. This result is at odds with the G09 
one so that it is worth investigating the cause. Compared 
to G09, there are three main differences. First, we have esti- 
mated gbar(ro) starting from the model parameters obtained 
by fitting the rotation curve of each galaxy rather than scal- 
ing the results in literature to the Burkert profile. While the 
scaling adopted by G09 typically works satisfactorily well 
(and, indeed, for the galaxies in the G09 sample we do not 
find a significant discrepancy), it nevertheless shifts some 
galaxies in the vertical direction thus possibly affecting the 
slope estimate. Moreover, our sample covers a smaller Mb 
range lacking low luminosity systems. Indeed, we find that 
the points missing in our sample, but present in the G09 
one severely departs from the extrapolation of our gb ar (ro) - 
Mb relation. Finally, we note that, while we only ask that 
a galaxy is well fitted by the Burkert model to include it 
in the sample, G09 applied some quality cuts before deriv- 
ing their results on the constancy of _w(ro)- Remarkably, if 
we apply their same cut, we find smaller slopes in marginal 
agreement with the G09 result. Such a circumstance makes 
us point at selection effects as a possible drawback of this 
kind of analysis. On the one hand, simulated samples may 
be used to fully mimic what is done with actual data and 
then investigate whether quality cuts may bias the inferred 
slope. On the other hand, one should increase the sample in 
order to directly check how the gt(ro ) - M B and g.(ro) - M„; r 
relations change depending on the selection cuts applied. 

As a final remark, we want to stress that it is worth re- 
considering the choice of the reference radius in this kind of 
analysis. As Fig.[T]shows, the estimated gi(R) are model in- 
dependent only if the reference radius R is carefully chosen. 
From this point of view, the scale radius ro is not an optimal 
choice since it probes the acceleration profile in the region 
where it is most model dependent if the halo concentration 
is too large. A preliminary study looking for a truly model 
independent quantity is therefore welcome before trying to 
understand whether the data argue in favour or against the 
existence of universal acceleration scale. 
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APPENDIX A: TABLES 

Below, we report relevant tables for readers interested in 
checking the analysis presented in this paper. Tab lefXT] gives 
the best fit, median and 68% confidence ranges on the fitted 
parameters (T^/T^j Cw, Vvir) for all the galaxies in the 
sample. Note that we do not directly fit the stellar M/L 
ratio, but rather its value scaled with respect to the fiducial 
one, T fid, adopted in the paper where the data are taken 
from. We also give constraints on the logarithm of the virial 
mass and the Burkert scale radius ro. 

Table [A"2l gives the data of interest for the fits involving 
the acceleration scales. We divide the table in two sections. 
First, we give the galaxies in the Good sample, then the 
ones discarded because of unreliable fit. We typeset in italic 
those systems belonging to the Gold subset and in italics 
those common to the G09 sample. We also report two dif- 
ferent quality flags. The first one refers to whether the rota- 
tion curve fit is reliable or not, while the second one checks 
whether the inferred ro is smaller or larger than Rmax, this 
latter being the last measured point of the rotation curve. 
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Table Al. Galaxy sample and constraints on model parameters. Columns are as follows : 1. galaxy id; 2. reference for the rotation curve 
data according to the following scheme: dBB02 = de Blok & Bosma (1992), S05 = Simon ct al. (2005), Sp08 = Spano ct al. (2008), 
THINGS = de Blok ct al. (2008), Oh = Oh ct al. (2011), Swll = Swaters et al. (2011); 3. reduced \ 2 = \ 2 /d.o.f. for the best fit model; 
4. best fit model parameters, 5. - 7. median and 68% confidence range for the fitted model parameters (Y*/Y f^, c v i r , V v i r ); 8., 9. median 
and 68% confidence range for the halo virial mass logMjji r and scale radius tq. Note that, for galaxies in the S05 sample, we have to set 
the stellar M/L ratio to its fiducial value since we have only data on the DM circular velocity. 
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Table A2. Quantities for the acceleration scales fits. Columns are as follows : 1. galaxy id; 2. absolute B band magnitude; 3. disc 
scalclcngth radius ; 4. RC quality flag (Qrc = for bad fit, Que = 1 f° r good fit); 5. Ro quality flag (Qo = for ro > Rmax, Qo = 1 
for ro ^ Rmax)', 6. - 9. median and 68% confidence range for loggi(ro) with i = H I ,d,bar, DM for gas, disc, baryons and dark matter 
and the accelerations in cm/s 2 . First part of the table refers to the galaxies in the Good sample, while the second one contains data for 
rejected systems. We typeset in italics systems belonging to the Gold sample and add a * symbol to identify those in common with G09. 
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